Low-T Phononic Thermal Conductivity in Superconductors with Line Nodes 
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The phonon contribution to the thermal conductivity at low temperature in superconductors 
with line nodes is calculated assuming that scattering by both nodal quasiparticles and the sample 
boundaries is significant. It is determined that, within the regime in which the quasiparticles are 
in the universal limit and the phonon attenuation is in the hydrodynamic limit, there exists a 
wide temperature range over which the phonon thermal conductivity varies as T 2 . This behaviour 
comes from the fact that transverse phonons propagating along certain directions do not interact 
with nodal quasiparticles and is thus found to be required by the symmetry of the crystal and the 
superconducting gap, independent of the model used for the electron-phonon interaction. The In- 
dependence of the phonon thermal conductivity occurs over a well-defined intermediate temperature 
range: at higher T the temperature-dependence is found to be linear while at lower T the usual 
T 3 (boundary-limited) behaviour is recovered. Results are compared to recent measurements of the 
thermal conductivity of T12201, and are shown to be consistent with the data. 

PACS numbers: 72.14.eb,74.72.-h 



Measurements of the thermal conductivity k in 
cuprates and other unconventional superconductors, 
made in the milliKelvin temperature range, have been 
valuable in studying T=0 properties of quasiparticles 
near the nodes of the gap. The measured heat current 
is carried by both electrons and phonons, and it is of- 
ten difficult to identify their respective contributions. To 
isolate the electronic contribution, the low-T data is typ- 
ically fit to a smooth curve that is extrapolated (eg. Ref. 
G]and references therein) to T=0. Assuming that the ex- 
trapolated value is dominated by electrons, the electronic 
part K e i/T(T = 0) can be obtained without knowing the 
temperature-dependence of the phononic contribution. 
This procedure has yielded values of K e i/T(T = 0) that 
are consistent with one another and with independent ex- 
perimental and theoretical result3222£££. However, if one 
wishes to obtain any temperature dependence of K e i/T, 
then some understanding of the phononic contribution 
is needed. This is particularly true when scattering of 
phonons by electrons is important, for in this case the 
thermal currents cannot easily be disentangled by ana- 
lyzing the magnetic field dependence of the data, ft is 
also in this case that the observed phonon heat current 
may be used to study the electrons. 

The importance of electron-phonon interactions in 
the low-T phononic thermal conductivity has been im- 
plicated in recent measurements^ on T^E^CuOe+a 
(T12201). The impurity bandwidth 7 is estimated to be 
10-40 K for this material, so the measured temperature 
range of below 1 K falls within the universal limit and 
K e i/T is expected to be independent of temperature. The 
data can be fit very well by an expression: 
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where a and (3 are independent of T. The phononic ther- 
mal conductivity K p t appears to be proportional to T 2 
(instead of T 3 , as expected for boundary-limited scat- 
tering) . Earlier measurementsiSii on overdoped BSCCO 



and recent studies^ of LSCO show similar behaviour. The 
authors of Ref. attributed the (3T term to phonons scat- 
tered by electrons, which is one obvious candidate for the 
relevant mechanism. 

Motivated by these data and by the general issues men- 
tioned above, 1 have calculated the thermal conductivity 
of phonons in the universal regime under the assumption 
that the significant current relaxation processes are scat- 
tering from nodal quasiparticles and scattering from sam- 
ple boundaries. The main result of this article is the fol- 
lowing: due to the presence of inactive phonon modes (i.e. 
modes for which the matrix element with nodal quasi- 
particles is zero) , there occurs a wide temperature range 
within the universal regime over which the phononic ther- 
mal conductivity n p h varies as T 2 . The width of this tem- 
perature range is such that the temperature may vary by 
a factor as much as A0/7, which can be between 10 and 
400 for optimally doped cuprates. The K p h oc T 2 be- 
haviour is expected to occur within a well-defined inter- 
mediate temperature range, with n p h oc T 3 predicted at 
lower temperatures (i.e. in the usual boundary-scattering 
limit) and n p h oc T predicted at higher temperatures. Be- 
low, I use experimental values for material-parameters to 
estimate the location of the temperature crossovers and 
find that they are consistent with the observation of T 2 
behaviour from below 100 mK up to several Kelvin in 
T12201. I have also compared the predicted magnitude 
of (3 to the value obtained from the T12201 data and 
found good agreement. 

Later in this article, I use a particular model for the 
electron-phonon interaction to calculate K p h when scat- 
tering by quasiparticles is important. The main qualita- 
tive result (the existence of a region in which n p h varies 
as T 2 ) does not depend on the model chosen, but follows 
from symmetry and other general considerations as I will 
first show. 

Symmetry arguments can be used to prove that the 
electron-phonon matrix element between certain acous- 
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tic phonon states and the quasiparticles at the gap nodes 
is zero (such phonon states are called 'inactive—). For 
example, it has been showni 3 - that in tetragonal crystals 
with a d x 2_ y 2 gap, a transverse phonon having its po- 
larization vector in the basal plane (the CuC>2 plane in 
cuprates) is inactive when the basal-plane component of 
its wavevector is parallel to the (110) direction. The in- 
teraction between inactive phonons and quasiparticles is 
obtained by expanding the matrix element, which varies 
with quasiparticle momentum on the scale of the Fermi 
wavevector kf, about the nodal points. In the universal 
limit, quasiparticles occupy a strip of the Fermi surface 
of length fcf (7/A0) centred on each node. The electron- 
phonon coupling for inactive phonons is thus smaller than 
that for active phonons by a factor 7/A0 and the corre- 
sponding phonon lifetime (proportional to the inverse- 
square of the matrix element) larger by Aq/j 2 . This can 
have a significant effect on low-T phonon thermal trans- 
port, which I now consider. 

Using a relaxation-time Boltzmann equation descrip- 
tion of the phonons, assuming that Mathiesson's law 
holds for the combined effects of scattering from quasi- 
particles and from the sample boundaries and ignoring 
any anisotropy in the speed of sound c s , we have 

n ph = C{k B T? I I dQ - m ^ i 

J dx> J 1 + (A/c s )tJ 

. (2) 

where T ph is the scattering rate of phonons by quasipar- 
ticles, A is the mean free path for boundary scattering, 
i.e. the length of the sample, C = AcJ 2 k B (2irh)~ 3 , n(x) 
is a Bose function and m{9, (ft) is the usual angular factor 
that gives heaviest weight to phonons moving parallel to 
the temperature gradient. This factor is not important 
unless the temperature gradient is applied perpendicular 
to the wavevectors of all inactive phonon states; it may 
be ignored in the present discussion. The first integral 
is over a dimensionless variable corresponding to phonon 
energy (only x between 1 and 10 contribute much to the 
integral), the second is over phonon direction. 

For sufficiently low phonon frequency, can be de- 
scribed in the hydrodynamic approximation, where it 
varies quadratically with phonon frequency and can be 
written as 

{Ncs)r^=x\T/T ff{<j>,e) (3) 

where To is a temperature scale and f{(ft, 9) is propor- 
tional to the electron-phonon matrix element. 

If the matrix clement is roughly isotropic, as expected 
for longitudinal modes, then the result for n p h will clearly 
be a function that crosses over at To from T 3 at low 
T (when boundary scattering, dominates) to linear-T at 
higher T (when quasiparticle scattering does). This does 
not adequately describe the T dependence of n p h in the 
latter case, since the main contribution comes from trans- 
verse modes for which f 2 (9, (ft) is highly anisotropic. 

When scattering by quasiparticles is dominant, the 
current is carried by the phonons that suffer the least at- 



tenuation by quasiparticles, i.e. inactive phonons. One 
may thus consider only modes for which states hav- 
ing particular directions of propagation will be inactive. 
Even for these modes, a typical direction of propagation 
will result in active attenuation (in the example given 
above, a transverse phonon polarized in the basal plane 
that propagates along [100] is active, so the lifetime of 
a phonon in this mode can be expected to decrease by 
a factor 7 2 /Aq as it is turned away from [110]). If the 
planar component (i.e. component perpendicular to the 
line node) of the wavevector of an inactive phonon makes 
an angle <fta with the x-axis then r~ h increases rapidly 
with — 00- The dependence on 9, measured along the 
node, is comparatively weak. 

There is a temperature region of width A0/7 over 
which the lifetime of active (inactive) phonons will be 
limited by scattering from quasiparticles (sample bound- 
aries). In other words, a region over which the second 
term in the denominator of Eq. 121 is negligible for (ft — (fto 
and large for typical (ft. Within this temperature region, 
the factor 1/[1 + (A/c s )t^ ] will be given by 



l + x 2 (T/T o ) 2 f 2 (9,0) xT\df(9,<P)/dcf>\ vr YUJ v ; 
so that 

k = C(k B T) 2 {k B T )r[4]([3]ir J dcos(9)n 2 (9, O ) (5) 

where I have combined the functions of angle into 
n 2 (9, (ft). The evolution of n p h with temperature is thus: 
T 3 at the lowest temperature, crossing over to T 2 at 
To (evaluated for an active phonon) then finally cross- 
ing over to linear- in-T at a temperature Tb(Ao/7). 

It should be noted that this result is valid in the hy- 
drodynamic limit of the phonon attenuation, in which the 
phonon-frequency dependence is quadratic. The origin of 
the T 2 -dependence is therefore unrelated to the usual T 2 
behaviour of this quantity that is seen in normal met- 
als, which occurs when the attenuation is in the quan- 
tum limit and has a corresponding linear-dependence on 
phonon frequencyi^i^. The T 2 regime is expected in the 
usual case for clean superconductors with line nodes, that 
is whenever inactive phonon states exist and 7/A0 << 1. 

Having established the existence of a T 2 range in k p u, 
I will now consider a particular model for the electron- 
phonon interaction in order to determine whether this T 2 
behaviour is likely to be observed under realistic condi- 
tions. This is done for the case relevant to cuprates so 
that a comparison may be made with the T12201 data. 

A form for the electron-phonon interaction in tight- 
binding square lattice materials such as cuprates is de- 
scribed in Ref. 0. The dependence of the matrix element 
on quasiparticle and phonon wavevector (k and q, re- 
spectively) are obtained for acoustic modes with the cou- 
pling strength characterized by the energy g = a\dt/dR\, 
the derivative of the hopping parameter t with respect 
to bond length in the effective single-band Hamiltonian 
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for the C11O2 plane (a is the lattice constant). Using 
this and, based on the discussion above, considering only 
transverse phonons polarized in the plane, the matrix el- 
ement is 



|M(k,q) 



(g^q) 

MNcI 



5(k,q) 



(6) 



where 



,g(k,q) = sin(0) 



A, 



• sin(20) + — cos(2(/)) 



M is the mass of the unit cell and N is the number of cop- 
per sites in the crystal. The angles are defined such that 
(6, <fi) = (7r/2,7r/4) describes a phonon travelling along 
(110), so (j) = tt/4. Also, g' = 4V2a\dt' /dR\ sm 2 (k f a), 
where kf is the Fermi wavevector at the node. The gap 
has been introduced using At/Ao = cos^a) — cos(k y a) 
to describe the k-dependence of the nearest-neighbour 
matrix element. Note that Eq. has the symmetry- 
imposed property mentioned in the previous section: it 
is zero for <f> = it /A and k x = k y . The first (second) term 
in Eq. 6 describes inactive (active) phonon states. 

The inverse-lifetime of phonons is equal to twice the 
magnitude of the imaginary part of the retarded phonon 
self-energy evaluated at the real phonon frequency cj q = 
c s |q|. I use the lowest approximation for the phonon 
self-energy, given by the 'bare-bubble' diagram with the 
electron-Green's functions taken in the universal limi1i^i& 
and the vertex given by Eq. [5] so that 
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Wq7V(k,q) 



I'll 



+ 7 2 )(jB 2 +7 2) 



(8) 



where I have used k B T, w q << 7 and h = 1. 

For active phonon states, the matrix element is inde- 
pendent of k and can be moved outside the summation. 
The contribution to the sum over k that comes from the 
vicinity of each node is given by 



E 



<i 2 



■n 



(£k +q + 7 2 )(^ + 7 2 ) 2v 2 k 



2Kf 



4 7 2 ^(aJ q /2 7 ) 



where 



S A {x) 



R(x) 



In 



1 + R(x) 



1 - R(x) 



(9) 



(10) 



R(x) = \/l — x 2 , no is the normal density of states, 
vikj w Aq, and o) q is the magnitude of the wavevector 
q = Vf<\ ■ ki + «2q ■ k2, with kj, and k2 unit vectors 
directed perpendicular, parallel to the Fermi surface at 
the node, respectively. 

In the hydrodynamic limilpi 3 " (Cj^/2^ — > 0), the factor 
in square brackets in Eq. OHis equal to w 2 . For large ujq/j 
(a limit validii only when c s q << 7 << v 2 q) this factor 
becomes 4^ ln(w q /7). 



For the inactive phonon state obtained by taking 9 = 
7r/2, = 7r/4 above so that g 2 (k, q) = A^/A 2 ,, the cor- 
responding result is 



E 



(A^/Ag)^ 7 2 



(^ 2 +q + 7 2 )(i? 2 +7 2 ) 



n 
2v 2 kf 



^|7 2 ^/K/27) 



with 



(11) 



(7) 5/(x) = 



2a; 2 7 2 

A2 

^0 



R(^) }S a (x) 



U:[^) + [l-R 

7 J \ x 

(12) 

where q~% is the k2 component. A phonon propagating 
with <p = 7r/4 has q 2 — for two of the nodes and = q 
for the other two. The latter give the main contribution 
to the sum over nodes since Si and Sa are monotonically 
decreasing with -Ek+q and Vf » v 2 in most cases. In the 
hydrodynamic limit the factor in square brackets in Eq. 
[TTlis (7 2 /A 2 ) ln(A /7)w 2 . This is the result noted above 
that, for phonons propagating in an inactive direction 
(<f> = tt/4:), the magnitude of r^ 1 is smaller than for 
typical directions by a factor 7 2 /Aq. 

To obtain n p h at temperatures well below T(A /^), 
I use the hydrodynamic limit for active phonons in 
Eq. [3 ignore (A/c s )t~^ for inactive phonons compared 
to unity and take the temperature gradient along the 
plane with m 2 (6», 0) = sin 2 (6>) cos 2 (^) and / 2 (6», </>) = 
cos 2 (2</))sin 2 (6 l ). The result is 



K p h = C(k B T) 3 J dxx 3 ( - ^ 



s{xt) (13) 



where r = T /To, s{x) — -^-[1 — (x 1 — x)atan(a;)], and 



(k B T y 



A 



( 9' 2 np \ 
TTc s v 2 k f \NMc 2 s )' 



(14) 



Eq. ^5] describes the smooth crossover, occuring at 
T ~ To, from T 3 to T 2 behaviour in K p h- In the limit 
T «T « Xb(A /7), Eq. US gives Eq. with the 
integral in the latter equal to ir/2. Note that both To 
and /3, the magnitude of the T 2 term in K p h, are uni- 
versal, i.e. independent of 7. Also, they do not depend 
on the detailed angle-dependence of the electron-phonon 
matrix element, but only on the fact that has a deep 
minimum at <f> — <f>o (see Ref. Il8j) . 

For a numerical estimate, I use g rs 6 eV obtained 
from observations of low-T downturns^ of K e i/T in opti- 
mally doped cuprates and assume that g' is of the same 
order— as g; the factor in parentheses of Eq. ^]is then 
of order unity and, introducing the Debye energy cd, 
k B To ~ y/ Aoeuia/A). For A equal to 0.3 mm (a value 
typical of the cross-sectional width of samples used in the 
measurements^^) , To is roughly 0.3 K and the coefficient 
(3 is 0.5 mW K -3 cm -1 . The former value indicates that 
it is plausible that the T 2 region may be observed within 
the sub-Kelvin T-range, and the latter value is within a 
factor of two of that measured for T12201. Also, using 9 
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7 = 3 meV and Ao = 40 meV, the linear-in-T behaviour 
would be expected to occur above 4 K (see Ref. l22|h 

The model gives similar predictions for other cuprate 
materials (doping dependence is considered below), 
whereas many data on cupratesA2222i indicate that K p h 
has a T-dependence that is closer to T 3 . However, the 
estimates used for the relevant parameters (especially g') 
are very rough, and the value of To obtained here is such 
that, within a T-range between 50 mK and 1 K or so, nei- 
ther T 2 nor T 3 behaviour can be ruled out. This work 
does provide a basis for understanding T 2 -dependence in 
cuprates when it occurs and should be considered when, 
for example, the magnetic-field dependence of k is ana- 
lyzed. A detailed comparison of this work with the low-T 
data is complicated by the observation of T 3 ~ 5 (where 
6 < 0.3) power laws arising because of a poorly under- 
stood boundary-scattering effeclsiS 3 - and because low-T 
downturns in n/T data occur due to heating losses^. 

I shall briefly discuss the doping dependence of these 
results. From Ref. 0, the T 2 coefficient (3 increases 
slightly with doping across the overdoped side of the 
phase diagram. Since To °c \J A /g' 2 , the results above 
are consistent with this dependence only if g' 2 decreases 
with overdoping, i.e. if it is roughly proportional to 
A . It is expected that g' , being derived from the effec- 



tive hopping Hamiltonian, should vary with doping and 
it would be interesting to see whether such a doping- 
dependence can be observed elsewhere. 

T 2 behaviour is less likely to be observed for under- 
doped materials. This is because Vijc a becomes very 
large (10 2 or so-) at extreme underdoping, so the hy- 
drodynamic approximation breaks down at temperatures 
much less than those for which the universal limit is still 
valid. (For overdoped samples V2/c s < 6 so the require- 
ment for the hydrodynamic limit, (v2/c s )ksT << 7, is 
not much more restrictive than that for the universal 
limit, k B T « 7.) If (c s /v 2 )j < T then (cJA)^ 1 
would not exceed one until temperatures much higher 
than To and, once it did, n p h would vary as T 3 /ln(T) 
rather than T 2 . This may explain why the T 2 depen- 
dence in Kph has only been seen for overdoped cuprates. 

In summary, I have found that the phononic thermal 
conductivity n p h varies as T 2 over a wide, experimentally 
accessible temperature range within the universal limit. 
This work may be useful in analyzing low-T thermal con- 
ductivity data, especially in studies aimed at determining 
temperature-dependent corrections to the universal value 
of K e i/T. I thank D. G. Hawthorn, Mike Sutherland and 
M. B. Walker for useful discussions. 
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